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V 

THEORY OF AN AEROPLANE ENCOUNTERING 
GUSTS, II. 

By EDWIN BIDWELL WILSON. 

i. This discussion is an immediate continuation of my previous 
treatment of the subject published in the First Annual Report of 
the National Advisory Committee for Aeronautics, Washington, 
X 9 X S> PP- 52—75 (Senate Document, 268, 64th Congress, 1st Session, 
reference to which will be by pages). The notations of that work 
will be continued without change except as hereafter noted. 

Periodic Longitudinal Gusts. 

2. That there is a certain degree of periodicity in gusts is obvious 
from casual observation, from the records of scientific observatories 
like Blue Hill, and from the familiar fact that all such phenomena 
in nature reveal a general tendency toward periodicity. Needless to 
say the periodicity is not mathematically exact in its regularity nor 
indefinite in continuance. 

The object, however, of an investigation of the effect of periodic 
gusts on an aeroplane can for practical purposes be no other than 
to reveal any exceptional effects that periodic, as compared with 
single, gusts may have upon the flight of the machine; and these 
exceptional effects will probably be indicated with sufficient practical 
completeness by an analysis built on the assumption of strict 
periodicity, long continued in operation — the phenomenon most to be 
feared being resonance. 

3. The longitudinal gusts are in, i°, head-on velocity u x ; 2°, 
vertical velocity w 1 ; 3 , rotary velocity q x . Very little is known as 
to the nature of rotary gusts (p. 65) and hence 3 may be left aside. 
It is not easy to see how vertical gusts can have any pronounced 
periodicity; the disturbance of the aeroplane's motion by vertical 
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WILSON— AEROPLANE ENCOUNTERING GUSTS. 213 

gusts is (p. 64) , except for very sharp gusts, essentially a convection 
of the machine with and by the gust ; for both these reasons 2° may 
be discarded. This leaves only i° — periodicity in the head-on 
gustiness — as likely to be of interest. 
The gust may be assumed in the form 

«! = / sin pt or «! = /*?*»". (1) 

The differential equations are (p. 59) 

f(D)u =— (0.12&D 3 + 1.160Z? 2 + 3.385Z? + o.9I7)m 1 , 

f(D)w= — Z> 2 (o.557£> + 2.458) w x , (2) 

f(D) $== — 0.02851!%!, 

with f(D) =D* + 8.49D 3 + 24.5D* + 3.385Z) + 0.917 

= (£> 2 + 8.359Z? + 23.37) (D 2 + 0.1308^ + 0.03924) . 

4. In the previous investigation it was found that the short- 
period heavily damped oscillation was not of much significance 
except in the case of a sharp up-gust (pp. 62-69), and that its 
significance in that case was not revealed in the major motion of the. 
machine but in the initial acceleration (or stress) upon it. It may 
therefore be expected that for periodic head-on gusts the short- 
period motion will be negligible in its effects. It is consequently 
desirable to carry out the numerical analysis in such a way as to 
separate, so far as may be, the short and long natural periods of 
the machine. 

Let us separate into partial fractions the operator 
_i 1 

f(P) ~ (£> 2 + 8.359D + 23.37) (K 2 + 0.1308Z) + 0.03924) ' 
or 

_i o.oi6£> + 0.089 — 0.01601D + 0.04263 

/CD) ~ £> 2 + 8.359Z) + 23.37 D* + 0.1 308P + 0.03924 • (3) 

The first fraction has to do with the short, the second with the long 
oscillation. The two operators are to be applied to certain ex- 
pressions derived from (1) by substitution in (2). 

5- If D = ip, the numerators of (3) have the respective magni- 
tudes 

(0.089 2 + q.016 2 ^ 2 ) 1 / 2 and (0.0426 2 + q.016 2 ^ 2 ) 1 / 2 . 
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For p = o, the second is about half the first ; for p — <x> , the two 
are equal ; the numerators therefore do not differ greatly in magni- 
tude for any value of p. 

The ratio of the denominators is 

T (.03924 -^) 2 +-i3o8 2 £ 2 l 1/2 
L (23-37 - P 2 Y + 8-329¥ 2 J ' 

and is very small when p is less than 1. For larger values of p, we 

have approximately 

i/i2 2 = i+2 3 .2//> 2 + 545/^. 

Hence the short oscillations may be neglected when p < 1 without 
introducing much error; but as p increases beyond the value 1, the 
importance of the short oscillation grows rapidly. 

6. Consider first the case p < 1, neglecting the short oscillation. 
The particular solutions for u, w, and 6, that is, I u , Ik, I e > are ob- 
tained from the imaginary parts of 

u — .01601 Ai + .04263 . .. „ „ 

-i = : ^r 1 . — - — — (.128* 3 ; + 1.16A 2 - 3.38^ 

/ -? + .1308^ + .03924 *■ r ^ v 6 6 ^ 

— .917V*", 

? ~ -01601^ + o 4 26 3 _ 2458^', (4) 

J - p 2 + .1308^ + .03924 v ° '^ ^ tsyj" > 

6 — .01601 Ai + .04263 . n .. . ( 

7 - - ^ + , 308^ + 03924 ( - - 028 5^)^ ( - 

To estimate the value of p corresponding to the maximum dis- 
turbance we may examine the amplitude of 9/ J, which is 

'd f (.o 4 263) 2 +(.oi6j>) 2 f 

amp -j - - 0285I 4 (.03924^+ (.1308^ j • (5) 

The calculation gives p 2 = 0.0394 or p — 0.1985. The value of the 
amplitude is then about 0.0095/ radians or 0.54/ degrees. If / 
should be 20 ft/sec, the forced oscillation would have an amplitude 
of about io°. 

7. As the use of />=* 0.1985 in calculating is somewhat more 
complicated than the use of /> = o.2, and as the change from 0.1985 
to 0.2 does not materially alter the amplitude of the forced oscilla- 
tion (and probably does not exceed the error of observations), we 
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may use p = o.2 in calculating the effect of a periodic gust of 
maximum resonance on the aeroplane. We shall first note that 
for /> = o.2 the ratio of the amplitudes of the two fractions in (3) 
is of the order 400 to 1, and the first fraction is therefore entirely 
negligible in determining the particular integrals. 
For the second fraction we have the complex value 

4.263 - .32* (4.275, - 4.3 ) 



-(££■-«*•)■ 



2.6i6i — .076 (2.617, 91. 6°) 
where the parentheses contain the polar coordinates of the complex 
numbers. The expressions into which this is multiplied to determine 
the coefficients of e*»* are for u/J, w/J, 0/J respectively 

— 0.922 — 0.6761= ( 1. 144, 216.24 ), 
0.0983 4-0.004561 = (.0984, 2.67 ), 

— 0.00571= (.0057, — 90 ). 
Hence the values of u/J, w/J, 0/J are 

«// = ( — .965 4- 1.651) (cos .2t-\-i sin .2t), 

w/J=( — .00918 — .1641) (cos. ,2t-\-i sin .2t), 

6/J=( — .00948 4- .000981) (cos .2t-\-i sin .2t), 

and /„ = 7(1.65 cos .2t — .965 sin .2t), 

/«, = /( — .164 cos 2t — .0092 sin .2t). 

I e = J(. 00098 cos .2t — .00948 sin .2t), 

I e ' = J( — .0019 cos ,2t — .0002 sin ,2t), 

J uo =i.65/, I wo =— .164J, I e0 = . 00098/, /'„„=— .0019/. 

8. On substituting these values to find the constants of integra- 
tion (p. 61), it is found that A and C, corresponding to the short 
oscillation in u, are negligible. Also B= — 1.65/, D = .J26J. 
Hence 

u=:Je-° Wit { — 1.65 cos .187*+. 726 sin .187O 

4-/(1.65 cos .2t — .965 sin .2t). 

In like manner (p. 62), A' and Q' are small and B' = .i76J, D' 
=—.051/. 
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w= /<?- 0654 <(.i76 cos .187* — .051 sin .187*) 

— /(.164 cos .2t-\- .009 sin .2t) — .oi2Je~*- 1Bt cos 2.43*. 

(The last term is added as a check on the initial condition w/ = o.) 
Finally (p. 62) y A" = .00007 J, B" = . 00104/, D" = .oiooJ, and 

=/*?-• 0654 '( — .00104 cos .187* + .0109 sin .187*) 

+ /(. 00098 cos .2t — .00948 sin .2t) + .oooo7/*r 4 - 18 ' cos 2.43*. 

9. Now to find the rise of the machine when the gust strikes it 
(p. 64). 

w + 115.56= Je-- 065it (.o$6 cos .187^+1.208 sin .187*) 

— /(.051 cos .2^+1.064 sin .2t). 
The cosine terms may be omitted. The integration then gives 
5=5.32/ cos .2f+o.44/ — Je-- 06Bit (2 sin .187^ + 5.76 cos .1870- 
A table of values of z may be computed as : 

t = o, 2, 4, 6, 8, 10, 12, 14, 

z/J = o, o, —.15, —.54, —1. 16, —1.90, —2.60, —2.97. 

This shows the rise or drop, according as / is negative or posi- 
tive, during the first quarter minute. The values of z now fall off, 
pass through o, and only become large as t nears 35. The natural 
oscillation is then becoming less effective relative to the forced 
oscillation which has a double amplitude of about 10.6/, or 202 ft. 
if / = 20 ft./sec. 

As the existence of a regular periodic gust for any long time is 
almost unbelievable, the only real interest in the calculation is in 
showing that during the first 15 seconds the effect of resonance fails 
to become so far established that the motion differs appreciably from 
that due to the simple head-on gust previously studied (p. 74). 

10. In the case of the machine constrained to remain horizontal 
during flight (by some automatic steering device), the corresponding 
equations (p. 69) are for u 1 = Je i v t 

u .I28pi + .598 



J -598 - £ 2 + 4-078^1 

w> = -557ft* 

J -598 - p 2 + 4.078/rf 



-. e %vt 



e %vt 
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As the natural motion is no longer periodic, there can hardly be 
any such thing as resonance, in the usual acceptation of that term. 
We can, however, ask what value of p will make w/J a maximum 
and hence induce the maximum oscillation in the vertical motion. 
To maximize 

p? I 

(.598 - p*)* + 4.078^ ° r 4-078 2 + (P - .59Z/PY 
take p 2 = 0.598 or /> = 0.774. The value of w/J is then 

w/J = — o. 1 36^*"', 

and the amplitude of w is 0.136/. The amplitude of the oscillation 
corresponding to the particular solution I w is 0.175/. 

Thus again it is seen that the steering device makes the motion 
far easier than when the machine is free (p. 70). There seems to 
be no need of carrying out the details of the integration. 

Note on Resonance. 

11. In defining, by implication, a state of resonance in the calcula- 
tions above, I have assumed that it was the angle which was to be 
maximized by the proper choice of the frequency p of the applied 
periodic force. It may be well to take up the theory of resonance in 
a little greater detail, for there are complications in the kind of 
system we have here to consider. 

A. G. Webster, in his "Dynamics of Particles, etc.," Teubner, 
1904, p. 175, gives general formulas for resonance and shows that 
if the damping coefficients are small and if the frequency of the im- 
pressed force nearly coincides with that of the natural oscillation, 
the amplitude of the forced vibration will be relatively large. 

This is not enough. For in the first place, the damping coeffi- 
cients in the case of the aeroplane can hardly be regarded as small 
(they sometimes exceed the frequencies) ; in the second place, we are 
not even certain that the motion of the system is wholly oscillatory 
(some of the roots may be real, and even positive if the machine has 
a certain amount of dynamical instability) ; and in the third place, 
under such conditions, the amplitude of the forced oscillation may be 
considerably greater when the frequency of the applied force is 
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materially different from that of the system (supposed oscillatory) 
than when the system and the force are nearly synchronous. 

12. The ordinary theory of simple resonance depends on the 
equation 

(D 2 + kD + n)x — J sin pt. 

The particular solution 

is the imaginary part of the expression 

Je ipt 



n — p 2 + kpi ' 

The amplitude of I x is the same as the modulus of the complex 
value x. The modulus of e ipt is i ; that of x is 

T I 

amp. i x [(w _£ 2)2 + /fe2 £2]i/2- 

13. To make the denominator a minimum we have merely to 
minimize 

O — q) 2 + k 2 q, q = p 2 >o. 

We find q = n — \k 2 , necessitating n> \k 2 . If, then, n > k 2 , the 
maximum amplitude of I x is 

max. amp. I x = 



where the positive or negative sign must be taken according as k is 
positive or negative. If n < |fe 2 , the maximum amplitude for I x 
occurs when p=o and is J/n. 

The amplitude is large when k or (n — ^fe 2 )* is small; it is very 
large when both conditions are satisfied. The largest possible value 
occurs when n=\k 2 and is \^2j/k 2 . In this case the applied force 
has an indefinitely small frequency where the natural oscillation has 
the frequency k/y/2. The theory of the system here considered is 
given by Webster (op. cit, p. 155). 

14. The case which corresponds to that in which we are in- 
terested is where the system starts from rest at the position of 
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equilibrium. The motion is then defined by the equation 



x = 



k(n- \k 2 ) 



( e -hkt cos Vw - \kH — cos Vn - J& 2 /) 



j / Vw — ?& 2 \ 

+ —, nr. ( , , - e _i4( sin Vw - |fc 2 * - sin <n - hkH ) . 

Under normal conditions this quantity remains tolerably small 
until the natural motion is nearly damped out or until that motion 
has time to increase greatly (k > o). Even if n = %k 2 -\-e 2 k 2 , the 
equation becomes 

2jt 2 J 

x = -p- (e -i *' cos |& — cos tkt) + -jj (te~ ikt sin \ht — sin e^/), 

and the conclusion still holds. 

For the simple system ordinarily treated for resonance the state- 
ment that the motion must be only slightly damped and the frequen- 
cies of the natural and forced vibrations must be reasonably near 
together, is therefore amply justified. The result holds even when 
n < \k 2 , in which case the maximum amplitude for l x (resonance) 
occurs when p = o and is J/n. 

15. The next simplest case is like that which arises in treating 
the constrained longitudinal motion (0=o) of the aeroplane (p. 69) : 

(D -\-a)u-\-bw= — au x — bw u a = .i28, b = — .162, 

cu-\-(D+ d)w——cu 1 — dw 1 , c = .557, 0^ = 3.95. 

The natural motion is given by 

A' =D 2 + (a + d)D + (ad — be) =0, 

and in this case by D 2 .+ 4.07&D + .598= o. Here the roots are 
both real, viz., — 3.93 and — 0.15. So far as the equation in D is 
concerned we have the case where k is large and n is small. The 
equations for the forced motion are 

A'« = — (aD +»)«! — bDw x , 

A'w = — (dD -\-n)w t — cu x . 

The question now arises: What is it that is to be a maximum? 
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For some purposes it might be the variables w or w — for example, 
the whole theory of gusts here given depends on the gust being 
small and producing. small effects, and if by an applied force, the 
values of u or w should become too large, the theory would become 
worthless. Again, if the question had to do with the strain on the 
machine, the derivatives du/dt and dw/dt would be the essential 
objects of interest, and should be maximized. Finally it might be 
the values x—fudt and z^fwdt — the actual displacements of the 
machine — which we desired to examine. Let us therefore consider 
several problems seriatim. 

16. Case i. — To maximize w with a head-on gust u t = e*"*. 

aip + n . .I2&ip 4- .598 

1! — fiipt = ! — pipt 

A' .598-^+4.098^ " 

The maximum value of 

.I28 2 £ 2 + .598 2 .I28 2 (£ 2 + 21.83) 



(.598 - p*Y + 4 .o 9 8 2 /> 2 p* + 15.59P 2 + -3576 
occurs when p 2 is o, that is, " resonance " occurs for p = o, the 
amplitude of the force and the oscillation being the same. 

Case 2. — To maximize w with a head gust. 

This was treated above (§ 10). The ratio .136 was found; the 
required value of p was .776. 

Case 5. — To maximize u with an up-gust w x . 

.i62pi 

u = e ipt 

.598 -p 2 + 4.og8pi • 

The condition is p =■ .776 as is Case 2 ; the ratio is .04. 
Case 4. — To maximize w with an up-gust. 

3.95jw +.598 .„, 

w = — — 6 P 

.598 - />» + 4.098/tf * 

The maximum value of 

3-95 2 £ 2 + -598 2 3-95 2 (^ 2 + .0228) 



(.598 - £ 2 ) 2 + 4.098*/* p* + i5.59£ 2 + -3576 
occurs when /> 2 =1.022 and p=.i$, and the amplitude ratio is 
about 1. 



WILSON— AEROPLANE ENCOUNTERING GUSTS. 221 

We note the very different values of p thus found, namely, o, 0.15, 
0.776, according to the choice of case. If in Case 1, we had taken 
/> = .I5, the amplitude ratio would have been about .7 instead of 1 ; 
if p = .776 had been assumed, the ratio would have been .37 instead 
of 1. 

Case 5. — If we desired to maximize 2 we should have had to 
treat 

1 3-95* + .598 1 P ipt 
i .598 -p 2 + 4.098^ e ' 

which would have given an infinite amplitude ratio for p=o. 

17. Now if we turn to the free machine and try to maximize 
fddt instead of 0, we have to maximize 

.04263 2 + .oi6 2 £ 2 



(.03924 2 - p 2 y + .1308 2 /) 2 

instead of (5, §6). The value of p 2 is about .0307 and of p about 
.175 instead of .2 as before. The amplitude ratio is then only 
slightly in excess (about 4 per cent.) of that previously found — in 
other words the numerical values are such that resonance for and 
for fddt, which is the preponderating term in the expression for z, 
occurs for considerably different values of p, but the effect is about 
the same. This may be regarded as validating our procedure (§6) 
in maximizing instead of f0dt. 

18. To sum up this discussion of resonance as applied to the 
aeroplane we may say that the frequencies which produce " reso- 
nance " depend largely upon the quantity in which the effect of reso- 
nance is to be sought and that the frequency which makes for a 
strong resonant effect in one quantity may make on another an effect 
much weaker than the maximum — or it may not. 

19. There remains to discuss the question whether the effect of 
resonance is practically serious, i. e., whether as in the case of the 
motion of the machine, above treated, the effect fails to make itself 
felt until after so long a time that the pilot would be entirely able to 
deal with it or the wind would really have in all probability ceased 
to be periodic with the period required. 

Now in order to insure that resonance is effective both of itself 
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and as against the natural motion, we should reasonably expect to 
require, i°, that the resonant frequency p be large (for if it be small 
the pilot will have ample time to take care of it), and that, 2°, it be 
reasonably different from any natural frequency which is only 
slightly damped (for in the latter case the initial conditions will 
probably be such as to cause the natural and forced effects to clash 
for a considerable interval of time). 

This problem in its generality is so complicated that I have as yet 
been unable to determine whether there may be practically serious 
effects due to resonance, but from the cases I have here treated, from 
the general considerations which I have advanced, with due regard 
to the restrictions on p which appear to be reasonable, and from cases 
which I have examined without mentioning them here, I should 
judge that resonance is not a practically serious matter in longi- 
tudinal motion, and that we may safely confine our attention to gusts 
of the form /(i — <r rt ). 

20. One type of resonance which deserves consideration is that 
of the damped harmonic gust Je~ nt sin pt. It would be conjectured 
that if — n± pi were nearly equal to a pair of roots of A = o, there 
might arise a considerable disturbance. It is not likely that a gUst 
of this type would exist in reality, but the commencement of any 
gust might resemble very closely the commencement of such a gust 
and if the effect of this type were very marked as compared to that 
of the types already considered, it would be necessary, for the sake 
of foreseeing the worst that could happen, to discuss this type. 

I have not time to take the matter up here. Moreover, I imagine 
that it would be found that the constants of integration turned out 
to have such values that the gust, though tuned in damping and in 
frequency to the natural motion of the machine, did not have very 
large effects except in cases where n and p were small enough to 
allow the pilot easily to correct for the disturbance. 

The damped periodic gust has been treated by Brodetsky, 1 who 
finds the amplitude of the particular solution is a maximum (for the 
machine I am dealing with) when t=i6 sec. and is then a tolerably 
large quantity, — but the pilot has a quarter of a minute in which to 
react to his environment. It is, however, by no means certain that 

1 Aeronautical Journal, London, 20, 1916, p. 154- 
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the pilot would have to react so quickly — the constants of integra- 
tion might turn out, as I have just suggested, such that the motion 
during the first quarter minute was not far different from that in 
the case of the simple gust. This was what was found to happen in 
the case of the periodic gust above treated (§9). The amplitude of 
the vertical motion so far as the particular solution was concerned 
turned out to be about 5.3/, but the constants of integration were 
such as to postpone the major effect of the particular solution until 
30 or 40 seconds had elapsed. If we have a damped harmonic gust 
and such a postponement were operative, the damping would become 
effective and the gust might turn out to have at no time an effect 
much in excess of the maximum effect of a single gust of the form 
/(i — e- rt ). 

Infinitely Sharp Gusts. 

21. In my previous paper I discussed gusts /(i — er rt ~) rising 
from zero to / with various degrees of sharpness depending on the 
value of r — the larger r, the sharper the gust. An infinitely sharp 
gust would be one for which r was indefinitely large. Such a gust 
would represent an absolute discontinuity in the velocity of the wind. 
This is impossible, though it represents a state of aerial motion 
which may be nearly approached. Moreover, the infinitely sharp 
gust could not strike the machine all over at once, and hence the 
theoretical effect of such a gust on the assumption that the machine 
is instantaneously immersed must differ from the actual effect upon 
a machine running into a discontinuity in the wind velocity. 

For this reason one may well limit his considerations to finite 
gusts with a value of r not greater than 5, say, as I did. Neverthe- 
less if the calculation of the effect of an infinitely sharp gust is 
simpler than for a finite gust and if the limiting motion derived for 
such a gust is not appreciably different from that for a sharp gust 
of reasonable sharpness, the discussion of the limiting case will be 
justified. 

22. Consider first the longitudinal motion and a head-on gust 
«i=/(i — e~ rt ), r enormously large. According to the symbolic 
method D = — r must be substituted to find the particular solution 
for e~ rt . As, however, A is of the fourth degree in D and all the 
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polynomials upon the right hand are of degree 3 or less, the result 
of the substitution is easy to find. 

For example, when u 1 = J(i — e~ rt ), 

I u /J=—i — e- rt (.i28/r), I uo =—J, 

/„//=— ■«-"(. 557/r), J«o = o, 

1,/J=> — ^02851/^) =0, I,o — o, 

I' e /J=e-'*(.028si/r2) = 0, /'*>=o. 

The equations of motion are 
M// = e- 4 - 18 *(.ooo9 cos 2.43* + . 0032 sin 2.43*) 

+ ^- 0654 '(.999i cos .187* — .3577 sin .187*) — 1 — e- rt (.128/r) , 
w/J = e- ilst (.io66 cos 2.43* — .0435 sin 2.43O 

_j_^-.0654*( — I0 66 cos .187^ + . 0352 sin .187O — ^" r '(.5S7A), 
1 000/7 = e~*- 1Bt ( — .0402 cos 2.43* — .0278 sin 2.43O 

+ e- 0864 *(.O4O2 cos .187?— .6683 sin .187O. 

The calculation of the constants of integration is much simplified. 
The terms e~ rt /r are retained because the stresses (forces) due to 
the gust are calculated from du/dt and dw/dt to which these terms 
make an initial contribution — there is an instantaneous initial stress. 
When t — o, 

du/dt = (.128 — .004 — .008 — .085 — .067)/= — .016/, 

dzv/dt= (.557 — .446 — .106 + .007 + .006)/= .018/. 

These are the initial accelerations and should vanish because the 
gust though infinitely sharp begins at zero. That they do not 
vanish is due to an accumulation of errors. 

23. Immediately after the initial instant, however, the first terms, 
viz., .128 and .557, being multiplied by e' rt vanish. The other terms, 
however, being multiplied by comparatively slow changing func- 
tions are not altered. Hence immediately after the first instant 
there are accelerations — .128/ and — -557/ along the x and z axes 
respectively. 



WILSON— AEROPLANE ENCOUNTERING GUSTS. 225 

To put it another way, there is a discontinuity in the stress at the 
initial instant — as might be expected. The amounts of the discon- 
tinuities are also just what might be expected, viz., X U J and Z U J. 
In like manner for an up-gust the initial discontinuities in accelera- 
tion are X W J and Z W J along the x and z axes. These results could 
have been foreseen from the differential equations themselves as 
well as from " common sense." The path in space is not materially 
different for an infinitely sharp gust from what it is for a reason- 
ably sharp gust. 

It may therefore be said that a tolerably good idea of what 
happens for sharp gusts may be had from the consideration of 
infinitely sharp gusts. 

24. It has just been stated that the conclusions concerning the 
initial accelerations may be foreseen from the differential equations. 
This may be proved as follows : We have 

( D — X u ) u — X w w — X q q — g$ = X u u 1 + X w w 1 + X q q u 

— Z u u -(- {D — Z w )w — (Z g -{-U)q = Z u u x + Z^w-^ + Zq<1>.> 

— M„u — M u w + (k B 2 D — M q ) q = M u u^ + M w w 1 + M q q ly 

De — q — o, 



(6) 



where the equations have been reduced to four involving only the 
first derivatives of the four variables u, w, q, 0, with the initial condi- 
tions u = w=q=0=o, by the device of choosing q = Dd as an 
independent variable so as to eliminate the second derivatives. 

These equations determine the first derivatives at the initial 
instant or at any instant in terms of the values of the variables at 
that instant, namely, 



Du =■ X„u + X w w + X q q + g0 + XuU^ + X u> w 1 + X q q u 
Dw = Z % u + Z w w + Z q q -\-Uq + Z tt u x + Z w w t + Z q q lf 

kjfDq = M u u -f- M w w + M g q + Af„« x + Mv/w^ -f M q q lt 
D6 = q. 
At the initial instant u, w, q, 6 vanish. 

PROC. AMER. PHIL. SOC, VOL. LVI, P, JUNE 21, 1917. 



(7) 
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With an infinitely sharp gust u lt w lt q t may be considered as not 
vanishing but as starting at finite values, J«, J w , J q . The derivatives 
are then at the initial instant 

UU = A-uJ u ~r* X-wJ w ~\~ XqJ q> 
U2V Z* uJ u ~\ £"wJ w ~\~ £ q-> q> 

(8) 

k B 2 Dq = MJ U +M W J W + Mqjq, 

Dd=o. 

The first two equations give the X and Z accelerations of the ma- 
chine which determine the stresses as the accelerations times the 
mass. 

We have, for numerical values, 

Du= — .128/* + .i62/ w -f-o/ g , if X q =o, 

Dw= — .S57J« — 3-95Jw + oJ q , if Z q =o, 

Z^Dq => o/„ + i .74/,,, — 1 50/,, if M u = o. 

The last equation determines the couple tending to break the ma- 
chine, by bending in the .ixsr-plane, on multiplication by the mass m. 

25. That which I have called an infinitely sharp gust is not an 
impulsive gust. The implusive gust is both infinitely sharp and 
infinitely intense, but endures for only an infinitesimal time. The 
effect of an impulsive gust is to produce instantaneous changes in 
u, w, q. Such an impulse, like the impulses of ordinary mechanics, 
puts an infinite strain on the machine for an infinitesimal time, and 
the only way to tell whether the machine will stand the strain is to 
take the yielding of the framework into account — it is a problem 
in elasticity. For the purpose of calculating the stresses produced 
by gusts on the machine I therefore prefer the sharp gust to the 
impulsive gust. 

For the purpose of treating the motion of the machine after the 
gust strikes it — the gust being now a sudden fierce squall in other- 
wise still air — we have merely to determine the constants of integra- 
tion from the initial condition u , w , q , and — o, where u , w , q 
are the impulsively generated velocities. These equations are 
(p. 61): 
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u = A + B, 

w =— 4.04A + 34-5C— .io58£> + .002587D, 

(9) 
0= — .132^4 — .0946C + .002478B + .005799Z), 

q = .703A + .205C — .001246.D + .000084D. 

Analytically the effect of the impulsive gust upon the equations 
for determining the constants of integration is merely to replace the 
initial values of the particular solutions I uo , I wo , I eo , I' m obtained 
on the hypothesis of finite gusts, by the respective values — m„, 
— w , o, — q . The effect of the disturbance may therefore be 
calculated at once from my equations (23), (24), (25), (26), as 
soon as the values u , w , q have been determined. 

26. In the calculation of w , zv , q the same doubt arises as in 
the theory of any very sharp gust, namely, the effect of the partial 
immersion of the machine. Is the effect of a blow traveling along a 
mechanism the same as that of the blow applied instantaneously at 
all points of the mechanism? The possibility of a difference be- 
tween the instantaneous immersion and the immersion distributed in 
time would arise only if, i°, the machine had time enough to change 
its orientation appreciably or, 2°, the acquired velocities were suffi- 
cient to change the relative wind and thus affect considerably the 
impulsive pressure. 

Even if we assume that no material difference in effect is to be 
expected, it is difficult to make the proper assumptions to lead to 
reasonably satisfactory values for «„, w , q for any actual machine 
whose characteristics are expressed in terms of the mechanical 
coefficients m, k B 2 , U, and the aerodynamical coefficients X u , X w , X q , 
Zu, Zw, Z q , M„, M w , M q . It is by no means certain that for a con- 
siderable aerial disturbance the finite instantaneous changes in u, w, 
q can be calculated from the equations (8) by replacing D by the 
sign A for the increment and taking /„, 7 W , J q as the intensities of 
the impulsive gusts ; for the nine coefficients X u , etc., vary with the 
intensity of the relative wind. 

It is for this reason that I have used finite gusts of various 
degrees of sharpness instead of impulsive gusts. Moreover, it is 



228 SYMPOSIUM ON AERONAUTICS. 

not certain but the finite gust represents more nearly actual condi- 
tions in the air when flying is at all possible. 

An article by Brodetsky, with an introduction by Bryan, has re- 
cently reached this country, 2 in which impulsive gusts are considered, 
relative to Bryan's skeleton aeroplane consisting of a forward main 
plane and rear tail plane. The discussion is both interesting and 
important as is everything to which Bryan, the great pioneer in this 
subject, sets his name, but it does not seem to help me, so far as I 
have yet been able to examine it, in regard to the effect of an im- 
pulsive gust upon a machine whose properties are actually de- 
termined in the wind tunnel. I have therefore decided to let stand 
the brief general considerations above. 

The Action of the Air Screw. 

27. In the work to this point, I have made for the discussion of 
gusts the same assumption concerning the action of the propeller 
that Hunsaker, Bairstow, and others have made for discussions of 
stability, namely, that under varying conditions the motor speeds up 
or slows down so as to deliver a constant thrust along the #-axis. 

It would be equally reasonable, from some points of view more 
reasonable, to assume that under changing conditions of relative air 
velocity a motor speeds up or slows down so as to deliver the same 
effective horsepower. We should then have the power P equal to 
the thrust H (taken positive) multiplied by the velocity — U : 

P = -HU= - (H + dH)(U+u), 
UdH + uH = o, 

u u 

<m=-n-=-p-. (10) 

This is an additional force which is directed along the X-axis if 
the propeller shaft is horizontal for the velocity of flight — U. If 
in the standard condition the shaft is not horizontal there would be 
components 

_ u u . 

— P — cos a, +P-™sma 

2 Aeronautical Journal, London, 20, 1916, 139-156. 
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along the x and z axes, a being the angle from the horizontal up to 
the direction of the shaft. Furthermore if the shaft did not pass 
through the center of gravity there would be a pitching moment 
— Phu/U 2 if h is the distance of the line of the shaft above the 
center of gravity. 

28. The equations for the natural longitudinal motion would 
then be 

(D-X u + ^Au-X w w-(X q D + g)8 = o, (11) 

the other two equations remaining unchanged, if we assume for 
simplicity that a = /i — o. The effect of the varying thrust is to 
change X u to X u — Pg/mU 2 . We have the value X u = — .128 for 
this machine. If the effective propeller horsepower were 87 for 
U= — 115. 5, the value Pg/mU 2 is 

Pg = 87X550X32 = q6 
mU 2 1800 X 13350 

The modification of the equations of motion on replacing X u 
= — .128 by X u = — .191 would make an appreciable, though not 
very serious change. 

The determinant A would become 

34D 4 + 290.8.D 3 + 850.9D 2 -f 165.1D -f 31.18 

= 34(D* + 8.553# 3 + 25.03^ + 4-^D + .917) 
as compared with 

34(1)* + 8.490O 3 + 24.50D 2 + 3.385^ + -917)- 
The rapidly damped oscillation would, as a first approximation, be 

— 4.276 ± 2.5961 instead of — 4.245 ± 2.5454. 
The first approximation for the small root would be 
— .097 ± .ijji instead of — .069 ± .i8ii. 

The damping would be more pronounced and the oscillation a trifle 
faster. 

29. It may be concluded that whether the screw is supposed to 
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deliver a constant thrust or a constant power is not very important 
to the theory either of stability or of gusts. It is not unlikely that 
the actual behavior of the screw lies within the limits set by these 
two assumptions or sufficiently near to one of the limits to validate 
the use of either hypothesis. 

The Aeronautical Journal, London, 20, 1916, p. 142, quotes 
Bairstow and Fage as giving the formula 

dH = — .onHdV, V in miles per hour, 

which is dH = — .ooj^HdV, V in feet per second. 

With [/= 1 15.5 numerically we would have for constant power 

dH= — .oo866HdV, V in feet per second; 

and, if I understand correctly the use of the signs -f- and — in the 
quotation, the results are in as good agreement as could be expected 
in view of the fact that I have no knowledge of the value of U 
for which the data quoted are given. (If the motor and screw were 
exactly designed to give a maximum efficiency at a standard speed 
U, we could not expect the efficiency to be the same at relative air 
speeds either higher or lower, and this would slightly influence the 
result.) 

Equations for Lateral Motion. 

30. The differential equations for the lateral motion of a machine 
in a gust may be written as (p. 54) : 

dv/dt + g<f> + Ur—Y v v—Y P p—Y r r=Y v v 1 + Y p p x + Y r r„ 

A/m. dp/dt — L v v — L P p — L r r=L v v 1 -\- L p p x -\- L r r lt (12) 

C/m. dr/dt — N v v — N p p — N r r =>N v v x + N p p t + N r r x , 

where the terms involving the small unknown product of inertia E 
have been neglected and gusts of the type v x , p lt r t have been 
allowed. 

The gust v x corresponds to a side wind. A change in the direc- 
tion of the wind by a small angle would produce such a gust even 
in absence of any change in the wind velocity. The gust p x is a 
rotary gust tending to produce a bank ; as a disturbance in the air 
it would correspond to a horizontal roller run into end-on (axially). 
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The gust r x corresponds to a column of air rotating about a vertical 
line. 

This last is a common type of aerial disturbance, easily observed 
on a warm day, often of very small diameter compared with the 
spread of the wings of an aeroplane, and accompanied by a strong 
rising current of air. Such a vertical vortex, if small, might strike 
one wing of the machine alone, and, due to the rising current, heel 
it over suddenly. It is, however, not this small local disturbance 
which we can consider by our methods here, but the larger and more 
gentle rotation in the air which might immerse the whole machine 
many times over and which produces a yawing motion in the 
machine rather than (primarily) a roll or bank. 

31. Place D = d/dt. Then the equations are 

(D—Y v )v+(g—Y p D)<f>+(U—Yr)r=Y v v 1 -\-Y P p 1 +Y r r 1 , 

—L v v+(k/D—Lp)b4—Lrr=L v v 1 +Lpp 1 +L r r 1 , (13) 

—N v v—N v D&+{.k c *D—Nr)r=N v v l +N P p 1 +N r r„ 

where k A 2 =A/m and k c 2 = C/m. The determinant whose vanish- 
ing determines the natural motion is 

D-Y v g-YpD 
A = - L v k/D 2 - LpD 
-N v - NpD 

Let the cofactors of A be 



U-Y r 

k c W-N r 



5n = 



612 = 



5l3 = 



622 = 



Sn = 



k/D 2 - LpD - L r 

-NpD k c 2 D-Nr 

kJD-N r lt| = _59 - 552? - 26 - 55 ' 

-N v 



= 2592P 3 + 23140D 2 + 8478D, 



k/D 2 - LpD 



-NpD 

- N P D k c 2 D - N r 
g-Y P D U-Y r 

D-Y v U-Yr 

- N v k c 2 D - N r 
-N v - NjJ) 

D-Y v g-YpD 



= - 32.84P 2 - 280.7D, 
= - 2270D - 868.8, 
= 70.6D 2 + 44.5D + 109.9, 
= 28.76, 
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^31 = 



g-YpD U-Yr 



k A 2 D* - LpD - L r 



= 4243D 2 + 36270D - 1776, 



U-Yr D-Y,\ 

Ltr X*k 



= 55.2.0 + 111.2, 



D-Y, g-YpD 
- L v k/D* - LpD 



= 36.7D 3 + 323.1I>* + 77 .88I> +27.15, 



where the numerical values are those arising from the data de- 
termined for the Curtiss Tractor (which is the machine under in- 
vestigation) by Dr. J. C. Hunsaker as given on page 78 of his paper, 
"Dynamical Stability of Aeroplanes," Smithsonian Misc. Collect., 
Washington, Vol. 62, No. 5, pp. 1-78, 1916, namely, 

Y v = — 0.248, Y p = o, Y r = o, 

£,,= + 0.844, L P = — 314, L r = + 55.2, 

N v = — 0.894, N p =o, N r = — 27.0, 

*/ = 36.7+. ft<?=70-6— , U=— 115.5, 0=32.17. 

The value of A is then (D — Y v )8 11 -\-g8 12 + C78 13 or 

A = 2592D 4 + 23780D 8 + 18000Z) 2 + 34610I) — 854. 

This result checks with Hunsaker's (loc. cit, p. 78) as well as prob- 
able. The equation A=o may be written as 

D* + 9.172D 3 + 6.943D 2 + 13.35D — 0.3295=0. 

32. From the last two terms, one root is indicated as D 
= 0.02468 ; and the correction can readily be found, giving 

D = 0.02436. 

There is another root near D = — 8.5, the exact value being 

D=— 8.542. 

The other factor of the biquadratic equation is 

D 2 + 0.6537^ + 1.583=0, 

of which the roots are 

D— — 0.3268 ± 1.2151. 
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The complementary functions for v, <j>, and r are therefore of the 
form 

v = C u e- 02486t + C 12 e- 8 "« + <r- 8268 '(C ]3 cos 1.215* 

+ C U sin 1.215O, 

<£= C 21 e-o<"* et + C 22 e- 8 - 642 ' + <?- 8268 <(C 23 cos 1.215* 

+ C 24 sin 1.215O, 

r= C 31 e 02iSit + C 82 e- 8 - 542 ' + e-™ 8t (C S3 cos 1.215* 

+ C 34 sin 1. 21 5O. 

The particular integrals for any gust may be represented as I v , 
I ^ I r , and their initial values as I vo , I^o, Ir , the derivative of /j, being 
I' 4 with the corresponding initial values 1'^. 

33. If as before (p. 59) we restrict the possible gusts to those of 
which the functional form is different from any of the four func- 
tions entering into the complementary functions, the particular solu- 
tions must, on substitution, annihilate the right-hand members of 
the differential equations, and the relations between the constants 
dj of integration may be determined from the two equations 

(D + 0.248)57 + 32.17.fr — ii5.5r=o, 

0.894^ + 04, + (70.6D + 27.0) r = o. 



Hence 

and 
Further 

and 
Finally 



•2724CU + 32.I7C 21 — 115-5^,1 =0, 
•894(7x1 + oC 21 + 28.72C 31 =0, 

C u =— 8.326C 21 , C 81 = .259iC 21 . 

— 8.294C 12 + 32.17Q2 — 1 i5-5C 3 2 = o, 
•894Q2 + oC 22 — 575-8^32 =0, 

C12 =■ 3-797C 2 2» C 32 = .005897^2. 



— .o788C 13 + i.2i5C 14 + 32.i7C 23 — ii5.5C 33 = 0, 

— i.2i5C 13 — .0788C 14 + 32.17Q4— 1 i5-5C 3 4 = o, 
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•894C 13 + 3-92C 33 + 8S74C M = o, 

.894C 14 — 85.74C,, + 3-92C 34 = o, 
and 

C 13 = i04iC 23 + 564.8C 24 , C 33 =— 6.37iC 23 + io.56C 24 , 

C 14 =— 564.8C 23 + io4iC 24 , C 34 =—io.56C 23 — 6.371 C 24 . 

The solutions therefore, so far as concerns the complementary 
function, are 

<t> = C 21 e<> 2i36t + C 22 e- 8s4it + e- 3268t (C 23 cos 1.215* 

+ C 24 sin 1.215O, 
^=— 8.326C 21 e 02 « 6 * + 3.797C 22 ^- 8 - 5 * 2 ' 

+ < ,- 3268 '[(io4iC 23 + 564.8C 24 ) cos 1<2ist 

+ (— 564.8C 23 + i04iC 24 ) sin 1.215*], 
r = 0.2571 C 21 e- 02436 ' + o.oo5897C 22 e- 8 " 2 ' 

+ «-» 1M *[(— 6.37iC 23 H-io.56C 24 ) cos 1.215* 

— ( 10. 56C 23 + 6.371) sin 1.215*]. 

34. These equations determine the relative magnitudes of the 
various sorts of natural motion. 

The first term is the slowly amplifying divergence, this machine 
being slightly unstable laterally. If a side gust is such as to induce 
a lateral velocity of — 8.326C 21 , it induces a bank of C 21 , an eighth 
as much in radians or seven times as much in degrees. It is 
therefore clear that only very small values of C 21 are admissible for 
safety. The second term, corresponding to the rapidly damped 
motion, shows such rapid damping that it can hardly be of impor- 
tance, except for possible strains on the mechanism, unless C 22 is 
so large that the whole work is inapplicable because of the failure 
of the motions to be small. 

The trigonometric terms show that the oscillation in v will be of 
great amplitude compared with that in <f>, the factor being about 
1200 when $ is in radians or 20 when <j> is in degrees; even the 
oscillation in r will be over 12 times as great as in <£. In other 
words, the machine may have a large oscillatory side-slip or angular 
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velocity of yaw without much bank, but for the divergent motion the 
bank is a serious matter for even moderate side-slip. 

35. The initial conditions <f>=:p = v=r=>o give 

O^t-21 -]- t-22 1 ^23 ~T * 4>o> 

0= .02436C 21 — 8.542C 22 — .3268C 23 + i.2i5C 24 + J'^ 

o= — 8.326C 21 + 3-797C 2 2 + i04iC 23 + 564.8^ + I vo , 

0= .2571 C 21 + .005897^2 — 6.37iC 23 + io.56C 24 +7 ro . 

These equations must be solved for the four constants C. 

C 21 =3 — .9839/ ^ — . 1 148/^0 + .ooo74O.T„ — .02797/ ro , 

C 22 =—. 000149/4,0 + .11701'^— .oooo3427 O0 — .on63/ ro , 

C 23 = — .01 595/4,0 — .0021 53/^0 — .ooo7o67 O0 + .O3967 ro , 

C 2i = .01468/4, + .ooi466/'4, — .ooo4537/ Vo — .07201/r,,. 

36. The equations from which the particular solutions are ob- 
tained are (since Y v =Np=Yr = o) : 

&v= (Z?8 U — A) v ± + LM .+ ( W» + NJ„)r v 

A<f,=DS 12 v 1 + L P 8 22 ^ + (L r 8 22 + .^2)^, (14) 

Ar = D8 13 v t + Ljdnp! + (L r 8 23 + NA 3 )r lt 
or 

Av= (— 640/?* — 9522D 2 — 34610D + 854K + (71340 

+•2732)^— (112560D 2 + no47ooD)r u 

A<t>= (—59-550 — 26.55)^ + (— 22150/? 2 — 13970.D 

— 345™)p l + (38950 2 + 97O0 + 3062)^, 

Ar= (— 32.84D — 28o.7)£> 2 z/ 1 + (—9030)/-! 

+ (— 992D 3 — 8724D 2 — 2103D +854)r 1( 
with 

A = 2592.D 4 -j- 23780D 3 + 18000D 2 + 34610/? — 854. 

Motion in Lateral Gusts. 

37. We shall take as before the type /(i — er rt ) for that of a 
single gust. 
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Case i. — Side-gust — sharp. v 1 =J(i — £"")• 

/,=/(— i + .oi473<r"), I.o=— .98527/, 
/$ = /( — .001028) <r 5 ', 7* = — .001028/, 

7'*=/(.oo5i4) < r 5t , /^ = .00514/, 

7 r = / ( .002706 ) e- st , I ro = .002706/, 

C 21 = — .000384/, C 22 = .0005364/, 

C 23 = .000809/, C 24 = .0002445/. 
The equations of motion are 

1000*//=— .384<?- 02 * 36 ' + .536e- 8 - 842 ' — i.028<r" 

+*- M68 *(.8oo, cos 1.215? + . 2445 sin I-2I50- 

This is all negligibly small. For the same reason certain terms may 
be neglected in v and r. 

*»//= .oo3<? 02436 ' + .oo2«r 8 " 2 ' — 1 + .oi473e- 5 * 

+ <?- 3268 '(.q8 cos 1.215? — .2022 sin 1.215O, 
ioor//=— .oi^- 02436 ' + .27ir 5 ' — r- 3268 '(.257 cos 1.215? 

+ 1.009 si n I.2I50- 
The effect of the sharp side-gust is to carry the machine side- 
ways with it, but not very powerfully at first — much of the air blows 
through the machine — the dominating term at first being 

t/= — .2/e - - 3268 ' sin 1.215?; 

after a few seconds the dominating term is v = — /, with the very 
slowly growing divergent term effective only after a considerable 
time. There is a slight yawing oscillation, but the extreme angle 
of yaw is only about 0.01/ radians or J/2 degrees — the angle being 
computed as 

iooiA// = f ioor/J-dt= .4(1 - e- mm ) + .054(1 -e~ u ) -.8316 

_j_ e --3268(( - 8 ; j I 6 cos 1.215?+ .0122 sin 1.215?). 

The actual sidewise velocity is compounded of v and the amount 
— 115.5^ due to the yaw. Hence 



y = J (v- U5.5t)dt. 
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For this calculation v and $ may be simplified to 

v = — / + /c 3268 '(cos 1.215* — .2 sin 1.215*), 

100^=— .378/ — 4/^- 02436 ' + /e- 3268 '(.832 cos 1.215O ; 
and 

3, = _. 5 67* + i8.5/(e 02436 ' — 1)— .146/ 

+ 7>- 3268 '(.i46 cos 1.215* + . 066 sin 1.215*). 

From this it will be seen that the oscillatory motion is, so far as 
concerns the lateral displacement, of very small amplitude. The 
first two terms which are progressive, are the ones which count. 
Moreover, the displacement is of the same sign as / although the 
side-slip v is of the opposite sign. This apparent contradiction is 
due to the change in orientation ^ — the machine moves away from 
the gust owing to the lateral excess wind-pressure, but turns into the 
gust owing to the moment of the pressures, and by virtue of the 
great forward velocity, this turning more than makes up, in the dis- 
placement, for the side-slipping. 

38. Case 2. — Side-gust — mild. t> 1 = /(i — c 2 '). 
/«,=/(— 1 + 1.0205c 2 '), ho = .0205/, 

7 ((l =/(.oo04043^- 2 '), 7^0= -0004043/, 

Z' 4 =/( — .0000809c 2 '), I '4,0= — .0000809/, 

7r = /(— .001514c 2 '), Ir =— .001514/, 

C 2 i = — .00033 1 /, C 22 = .00000738/, C 23 = — .0000807/, 
C 2i = . 0001055/. 

It is again seen that there is practically no rolling motion pro- 
duced by the side-gust. For v and r, 

v/J=.oo27e 02i3et — 1 + 1.0205c 2 ' 

_l_<,-.3268«( — 0244 cos 1.215* + . 1554 sin 1.215*), 

ioot/J = — .oo85^ 02436 ' — .1514c 2 ' 

+ c 3268 '(. 1628 cos 1.215* + -.0672 si" 1-215*). 

(The check v = o, r = o, when * = o, shows that the accuracy has 
been reduced so that the third place is not sure,) The effects of 
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the gust are qualitatively as before. The oscillatory motion is not 
pronounced ; the ultimate side-slip velocity is — / ; the ultimate dis- 
placement has the same sign as / because the divergent term in 
v — n5-5i> is positive. 

39. Case 3. — Side-gust — oscillatory. When one examines the 
records made or making at such an observatory as Blue Hill for 
gustiness in the air, no phenomenon is perhaps more striking than 
the reasonably periodic side-switching of a reasonably steady wind. 
A south wind, for example, may whip back and forth between 
S.S.E. and S.S.W. for hours at a stretch, as Prof. Alexander 
McAdie has been kind enough to show me on some of his records. 
In the absence of rotary motion, concerning which I am unable to 
find satisfactory data, the simplest way to figure this change in direc- 
tion is as a periodic side-gust. A machine going south in such a 
wind would experience an alternating side-gust. (The oscillations 
in the head-on velocity of the wind would be relatively very small 
except for actual changes in head-on velocity superimposed upon the 
changes in direction.) It is therefore especially interesting to dis- 
cuss a periodic side-gust — this being the only periodic gust of which 
we can reasonably be said to know anything at all definite. 

Let v 1 = Je ipt . We may assume, from our work above that the 
rolling motion will be small and that the side-slip velocity v will not 
be of as much importance in determining the path as the angle \p 
coupled with the large forward velocity. The complex value of r is 

(280.7 ~ 32.HPi)P 2 Je ipl 

_ 2592^ - i8ooo£ 2 - 854 + ^(34610^ - 23780^) ' 

If at any one place the period of the complete oscillation is 2w/n 
with the wind velocity V, the distance traveled by the wind during 
the time of an oscillation in direction is 2irV/n, and this is the dis- 
tance between the nodes of the motion. The time required for this 
machine (£/=■ — 115. 5) to pass over the distance 2-n-V/n is 
2-rrV/ii$.$n. The periodicity of the gust as it appears to the 
operator of the machine will therefore correspond to the value p 
= ii5.5w/F. For instance, if F = 20 and the time of an oscilla- 
tion at one spot were 10 sees, so that 71=0.63, the value of p would 
be about p = ^.6, and the oscillations would appear to the pilot as 
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taking place about every if seconds. A slower oscillation, i. e., a 
longer periodic time, would diminish n and p, — an oscillation at one 
spot every half minute corresponds to a value p = i.2 on the basis 
of the assumptions made above. 

In considering the values of p which make the amplitude of r 
large, the only hope is to make the term 34610/) — 2^y8op s tolerably 
small. This means p 2 must be about 1.5. For this value, the 
modulus of r is about .03/ and the modulus of the yawing oscillation 
corresponding will be about .025/. If a wind of 20 ft./sec. is 
whipping through an angle of 45 °, the side-gust will be only of 
about 7 ft./sec. semi-amplitude and the angle of yaw will be in the 
neighborhood of .175 radians or io°. There is nothing to indicate 
that this would be fatal, though it would surely be a nuisance. 

Owing to the fact that the coefficients of i in both numerator and 
denominator are relatively small, the angular velocity I r would be 
about in phase with the gust v ls and hence the angle 1^ would be 
about quartered in phase. If there were periodically an angle of 
io° or 12 between the direction of flight and the relative wind, we 
should find that we were getting into a region where considerable 
rolling and pitching might be induced — for as Hunsaker has shown 
(loc. cit., p. 62) the lateral and longitudinal motions are not strictly 
independent; but as the machine makes the major part of the rela- 
tive wind, the directions of flight and of the relative wind never 
differ greatly — only some three degrees at most in the case under 
consideration. 

It seems hardly necessary at this time to go into the calculation 
of the actual motion; enough has perhaps been accomplished in 
showing that the oscillation of the direction of the wind induces at 
most a moderate yawing of the machine. The semi-amplitude of 
115.5^ would be, if 7 = 7 ft./sec, about 20 ft. ; the center of gravity 
of the machine would sway back and forth across the line of flight 
with a total amplitude of 40 ft., until the divergent term became 
effective. 

40. Case 4. — Rolling gust. p 1 = J(i — e~ rt ). If there were no 
interaction between v, p, r, the effect on rolling of a rolling gust 
would be figured from the equation 
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36.7DP + 314P = - 3i4->Xi - e~ rt ), 
pi J = - 8.055<r 8 - 055 ' f e 8 - 055, (i - er")dt, 
p/J= . I+e -s,, t+ 8 -°55 ,, 8.055 



5.055 - r 8.055 - f 

This means that for any ordinary sharp gust p rapidly acquires the 
value — /, and <f> the value — Jt (radians). It must therefore be 
expected that unless / is very small indeed, the motion will be much 
disturbed. There will be developed a component of the weight in- 
ducing side-slipping, and yawing will rapidly develop — the machine 
apparently goes off on a spiral dive. 

We may make the calculations in detail when r= 1. Here 

V/=— 3.14— .114^"*, ho=— 3- 2 5^ 

7^/7=40.4— • 1. \er\ /«o = 39-3^ 

/y/=i.i<r<, /' #0 =i.l/, 

I r /J= 10.58— .234<r', I ro — 10.35/. 

C 2l =— 39.1/, C 22 = .oo27/, C 23 =— .219/, C 24 =— .163/. 
The equations of motion become 
<£//=— 39.ir 02436 ' + .oo3e- 8 " 2 ' + 40.4— i.i<r* 

_|_<?- 3268t ( — .219 cos 1.215* — .163 sin 1.215*), 
sy7=324£- 02436 ' + .oiO3*r 8 - M2t — 3.14 — .1140"' 

_j_£-.326«( — 220 cos 1.215* + 49.2 s i n i- 2I S0» 
r/J=— ioe 02i36t + 10.58 — .234T-' 

+ e- 326S '(— .33 cos 1.215*4-3.35 sin 1.215*). 
In the equation for </> the effective terms are 

4>/J=— 39(r 02436t — 1) =— f(nearly), 

and there is a steady divergence in <j> to the approximate amount 
— /* as foreseen. The side- ways velocity v develops more slowly, 
perhaps, but after one second amounts to something like 300/. It 
is clear that / must be very small or the motion becomes disastrous. 
It would be of especial interest to know what sorts of magnitudes 
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for / are likely to arise in flight under normal conditions. In so far 
as experience shows that machines are not liable to roll and side- 
slip, it is pretty good evidence that aerial rotary motion with axis 
parallel to the earth is rare and small. 

41. Case 5. — Yawing gust. ^=7(1 — e-*). 

I v /J= 25.670-*, / ro = 25.67/, 

V /= — -°792 + .154^, 1 40= -075^ 

I\/J=— .154*-', 7' 4 o=— .154/, 

Ir/J = — 1 — .I235<r«, ir = — .1.12/, 

C 2 i= — .006/, C 22 = — .006/, C 23 = — .0634/, C 2 4 = .0702/. 
In this case the motion is 
4>/J=— .oo6e 02iWt — .oo6<r 8M2t — .0792 + .154*"' 

_j_ e -.3268f(- — .0634 cos 1.215* + . 0701 sin 1.215*), 
v/J = + .o5e- 02436 ' — .023<r 8 - 542 ' + 25.67*?-' 

_|_ e -.3268*( — 26.35 cos 1.215*+ 108.9 sin i- 2 i5*)> 
r// = — 1 — .1235^' + <?- 3268 '(i. 145 cos 1.215*+. 222 sin 1.215*). 

For moderate values of /, there is nothing serious indicated. The 
coefficients of the divergent terms are small. There cannot be much 
roll. The most noteworthy phenomenon is the large amount of 
side-slip which is fairly rapidly damped out. 

42. This leaves the rolling gust as the only dangerous type of 
lateral gust. 

The infinitely sharp side-gust would produce an initial accelera- 
tion Y V J. 

Constrained Aeroplanes. 

43. Suppose now that by some automatic steering device the aero- 
plane were constrained to remain pointing in the same direction, 
i. e., so that r = o identically. The equations of motion become 

(D - Y v ) v + {g - Y P D) 4, = Y v v, + Y p p t + Y r r u 

— L v v+ (k A 2 D — L P )D<j> = L v v 1 + L p p 1 + L r r 1 , (15) 
— N v v — N p D<j> = N v v x + N p p 1 + N r r t + F, 

where Fm is the moment necessary to maintain the constraint. The 

last equation may be regarded as determining F. 

PROC. AMER. PHIL. SOC, VOL. LVI, Q, JUNE 20, I917. 
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The natural motion of the constrained machine is found from the 
determinant 

A' = 833 = 36.7£ 3 + 323-iI> 2 + 77-88D + 27.15 = 0. 

This is a cubic equation which has no positive root. 

The negative root is — 8.54. The quadratic factor remaining 
after division by D + 8.54 is 

3 6.7£ 2 + 8.746D + 3- 18=0, 

of which the roots are 

D = — 0.1 19 ± 0.2691. 

The real part is negative and hence the motion is dynamically stable. 
The introduction of the automatic device has removed the in- 
stability in the lateral motion. As compared with the complex roots 
in the free motion, these roots indicate a much slower period and a 
considerably smaller damping. 

44. On the other hand suppose that the constraint had been such 
as to keep the machine level, i. e., <j> — o identically. The equations 
would have been 

(D— Y v )v+(U — Y r )r= Y v v t + Y p p 1 + Y r r x , 

— L v v — L r r=L v v 1 -\-L p p 1 -\-L r r 1 -{-F > (16) 

— N v v + {k c 2 D—Nr)r^N v v t + N p p 1 + N^. 

The natural motion would have been determined by 

A" = S 22 = 70.6D 2 + 44. 5D + 109.9 = o. 

The roots are 

P = — 0.31 5 ±0.2371. 

The machine is again stable. 

45. It follows that at high speed this Curtiss Tractor, which is 
laterally unstable when free, becomes quite stable when constrained 
either to remain on its course or to fly on even keel. 

If stabilizers against rolling and turning were provided, the 
motion would reduce to 

(D — Y v )v = Y vVl + Y,p x + Y r r lt (17) 

and would be stable, D— F„=> — 0.248. 
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46. It would be a relatively easy matter to discuss the effect of 
gusts of various types on the aeroplane constrained in various ways ; 
two equations are much easier to handle than three. Until some 
definite problem is proposed as important, until some particular con- 
straining device is indicated as likely to be adopted, it may be as well 
not to go into the calculations, which are quite straightforward. 

That a constraint against rolling might be worth while, and 
would indeed be very valuable if rolling gusts were a common thing, 
is suggested by the work done on the free machine (§42) where 
gustiness was seen not to be very serious except for the rolling gust. 

Discussion of Method. 

47. I pointed out in my earlier paper that there were several 
outs about my method of treating gusts. First the gusts must be 
small. If they are not tolerably small, flying would be too difficult — 
so that assumption is not wholly unjustifiable. Second, the calcula- 
tions for determining the individual equations of motion and for 
determining formulas for the constants of integration are very 
tedious. Third, the numbers are of such various magnitudes that 
the arithmetical operations which must be carried out cut down the 
accuracy of the work a good deal and indeed, unless great care is 
taken, will lead to illusory or incorrect results. This does not ap- 
pear to be due to any very rapid variation of the true results cal- 
culated from varying data, but to the mode of computing. 

To offset these inconveniences we have the satisfactory result 
that once the preliminary calculations are made, many and varied 
types of gusts may easily be treated, and the further valuable result 
that the actual motion for each case is known so that not only the 
initial motion is determined, but the whole extent of the motion. 
This last is necessary for any just appreciation of the effects of 
periodic gusts and resonance, as has been shown. 

48. For another method of treating gusts reference may be made 
to a recent paper by Brodetsky and Bryan, " The Longitudinal Initial 
Motion and Forced Oscillations of a Disturbed Aeroplane," Aero- 
nautical Journal, London, 20, 1916, 139-156, which has already 
been cited in the text. 
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Much may be said for their method. of expansion in series — for 
some problems the work is decidedly simpler than with my method. 
It has been my experience, however, that the application of series 
to the motion of any aeroplane has its own difficulties and com- 
plicated calculations when the motion is to be followed for any 
reasonable length of time and especially if the machine is defined, 
as I have always preferred to regard it as defined, by the actual 
coefficients determined by wind tunnel experiments rather than as 
Bryan's skeleton plane consisting of a main front plane plus tail 
plane, — even though the results obtained from such a skeleton may 
be extended to more complicated machines by Bryan's invariant 
method (see his " Stability in Aviation," Chap. VI.). 

49. The question therefore arises whether there may not be 
some way of abridging the calculations leading to the actual motion 
of the machine. Since finishing my work above, I have received the 
Proceedings of the London Mathematical Society, 15, 1917, Pt. 6, 
in which there is an article on " Normal Coordinates in Dynamical 
Systems," by T. J. I'A. Bromwich in which he develops a method 
of treating the motions of dynamical systems by means of the theory 
of functions of a complex variable. I wish, in closing, to describe 
the application of Bromwich's work to the problem in hand. 

We have to solve for the longitudinal motion equations of the 
type 

(D — X u )u — X w w— (X q D+g)6 = P 1 e* t , 

Z u u+(D — Z w )w—(Z q +U)DB = P 2 e» t (18) 

— Muu — M w w+(k 2 D 2 — M q D)0 = P i e* t , 

B 

where p is a real or complex number, the values we have used being 
o, — r, ± pi. We substitute 



1 f 

U = : I 1 



e A, #X, 



w = ^-.\ e M nd\ (19) 

2in J 



--■(■ 



e"?d\, 
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where the integrals are loop integrals in the complex plane and 
I, t), f are any functions of A. The results are 

^ jf [(X - X U )S - X wV - (X t \ + gMe M d\ = P^\ 

Tlfi- Z*i + (X - Z a ) v - (Z, + C7)Xf]e x< JX = P 2 e* ( , (20) 

Tlfl- Mu$ - M aV + (k B \ 2 - M g \)fte"d\ = P^\ 
We next set 

(X - X u )t - X aV - (X t \ + g)f = Pj/CX - /i), 
- Z u t + (X - Z„)i» - (Z, + C/)Xf = P 2 /(X - /.), (21) 

- M M ? - Af„i, + (fe B 2 X 2 - M s X)f = P 3 /(X - /i), 
and solve for |, 77, f, finding 

Pi^ii + Pidn + P3531 



t = 



A(X - M ) 

PA 2 + P 2 5 22 -\- Pzfa: 
A(X - ix) 

PiSis + P 2 5 2 3 + P353; 



(22) 



A(X - M ) 
A = 34(X 4 + 8.49X 3 + 24.5X 2 + 3.38SX + .917) . 

Bromwich shows that, if with these values of £, y, £ we take the 
loop integrals (19) around a very large circle, the results for u, w, 
9 will be the solutions for the motion disturbed from rest at the 
position of equilibrium by the impressed forces P. As he points 
out, this integration is equivalent to the sum of the integrals around 
infinitesimal circles about X = /t and about each of the roots A of 
A=°o, that is, the integral is equal to the sum of the residues of 
ge**, t)e xt , t,e xt . There is no need to calculate any constants of in- 
tegration. Moreover any of the quantities u, w, $ can be obtained 
without the others. The numerators in g, r/, £ are already calculated 
in (20 a, b, c) of p. 59. 
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We have, for example, for a head gust u u 

■ i28\ 3 + i-i6\ 2 + 3-385X + -9i7 
* (\-M)(X + 4-i8±2. 4 3i)(X + .o654±.i87i) Ml ' {23) 

where the double sign stands for two factors, and u 1 ^=J{\ — e~*), 
to take a particular case. The residues at each point are merely 
the values of the fraction when one of the factors, the one which 
vanishes at that point, is thrown out of the denominator. In the 
first case for i = £°* we have as residue of £e xt = $: 
at \ = fi=o, 

W 

(+ 4.18 ± 2. 43 i)(.o654 ± .187^) 

at A = — 4.18 — 2.431, 

-I28X 3 + i.i6X 2 + 3-385X + -917 

(— 4.18 + 2.43i)(- 4.86i)(- 4.12 - 2.43; ± .187*) ' 

at k= — ; 4.18 + 2.431, the conjugate imaginary expression. And 
so on. To treat er % we should have: 
at A.=/u= — i, 

_ - .128 + 1. 16 - 3.385 + .917 

(3.18 ±2.43*) (.9346 ±.187*) ; 
and so on. 

As the calculation with imaginaries involving squares, cubes, 
products, and quotients is by no means simple, it is clear that to get 
the solution for u will be reasonably hard work — much harder than 
to find the particular solutions which for the simple gust involved 
only, real numbers. It may be admitted that to work any one gust 
the labor will probably be much less than by my method of determin- 
ing formulas for the constants of integration in terms of the initial 
values of the particular integrals. But as far as I can see, Brom- 
wich's method is of no particular advantage if we desire to calculate 
the effects of a large number of gusts /(i — e~ rt ) of various degrees 
of sharpness both head-on, up, and rotary. When we came to cal- 
culate a periodic gust we found that we were involved in powers and 
products' and quotients of complex numbers, and it is probable that 
the work we did in finding the particular integrals was comparable 
with that required for the present analysis. 
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Summary. 

In continuation of my previous work in gusts as affecting the 
Curtiss Tractor JN2, I have discussed : 

1. Periodic Longitudinal Gusts. — It was found that, even in the 
case of best resonance with the slow natural oscillation, the motion 
was not much different from that produced by a simple head-on 
gust until after a considerable time (over 14 sec.) had elapsed. The 
amplitude of the forced oscillation (in up and down motion) which 
ultimately became effective was about 5 times the amplitude of the 
gust. This was not regarded as serious because true periodicity can 
rarely be maintained in a head gust and because no pilot would 
wait to let its effect reach such a magnitude. Periodic up gusts and 
rotary gusts were considered as not likely to arise. 

2. General Theory of Resonance. — It was shown that for aero- 
plane problems resonance meant different things for different prob- 
lems. It was inferred that resonance was unlikely to be particularly 
serious because in all probability its effect would either be small 
or would take so long .to become established that the pilot would 
check it. 

3. Infinitely Sharp Gusts. — It was seen that the shock to a ma- 
chine was mX u J and mZJ for a head gust, and mX w J and mZ w J 
for an up gust. The serious case is mZ w J, the vertical shock in an 
up gust which was about $J/g times the weight, more than twice that 
found for the sharpest gust previously treated. It would be still 
more serious in a machine where Z w was greater than in the JN2. 
The Moral: Keep Z w small, clashes with Hunsaker's conclusion 3 
that lateral stability is incompatible with high wing loading — but 
such an antithesis is common. 4 Reference was made to impulsive 
gusts. 

4. The Effect of the Propeller. — The assumption that a constant 
power instead of a constant thrust was delivered did not very ma- 
terially alter conditions of flight. 

5. Lateral Gusts. — The general equations were set up and 
integrated. 

3 " Dynamical Stability of Aeroplanes," Washington, Smithsonian Misc. 
Collect, 62, 1916, p. 77. 

* " The production of a laterally stable aeroplane is attentant with many 
compromises," Hunsaker, p. 74. 
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(a) Single side-gusts were shown to produce modern side- 
slipping, insignificant roll, and moderate yaw. It was seen that the 
yaw was into the relative wind so that the displacement of the ma- 
chine in space was toward the gust despite the side-slipping. 

(b) Oscillatory side-gusts were shown to be a common condi- 
tion of flight, to produce moderate side-slipping and yawing, but 
insignificant rolling. The path of the center of gravity proved 
to be sinusoidal, so far as the forced oscillation was concerned, and 
of amplitude about 2 or 3 times the amplitude of the gust. 

(c) Yawing gusts were found to induce a good deal of side- 
slipping, but did not appear to be serious. The roll was very small. 

(d) Rolling gusts were seen to put the machine into a spiral 
dive, and thus to cause a real danger unless the motion were checked 
promptly by the pilot. 

6. Constrained Machines. — A device to keep the aeroplane on its 
course or to prevent rolling made the previously unstable machine 
stable. Such a device might be important to reduce the liability 
to the spiral dive in rolling gusts provided such gusts were common 
phenomena in flying weather. 

7. Other Methods of Treatment. — The Bryan-Bordetsky method 
of initial motions and Bromwich's new method of finding the solu- 
tion for a disturbed state without calculating the constants of integra- 
tion were briefly compared with my system of analysis. 

Massachusetts Institute of Technology, 

Cambridge, Mass. 



